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An  analytic  solution  to  the  optimal  low  thrust  transfer  problem  is  developed  using  e 
first-order  perturbation  approach  for  planar  orbit  transfers.  This  solution  is  verified  by 
comparison  with  the  solution  obtained  by  integrating  the  equations  of  motion.  For  low  vehicle 
accelerations  the  analytic  solution  yields  results  quite  close  to  those  obtained  by  the  integrated 
solution.  At  higher  acceleration  levels,  however,  the  analytic  solution  diverges  from  the 
integrated  solution.  For  very  low  acceleration  levels  the  analytic  solution  approachs  the 
limiting  cobb  of  the  infinitesimally  low  thrust  spiral  solution. 


OPTIMAL  LOW  THRUST  ORBIT  TRANSFER 
USING  A  FIRST-ORDER  PERTURBATION  MODEL 


I.  Introduction 


In  performing  any  orbital  transfer  it  is  essential  to  determine  the  opttmum  thrust  vector 
profile.  This  problem  is  easily  solved  for  two  extreme  cases:  the  continuous  thrust  spiral 
transfer  far  a  vehicle  having  infinitesimally  small  acceleration,  and  the  two  burn  impulsive 
transfer  for  very  htyt  vehicle  acceleration.  Between  these  two  extremes,  the  optimal  transfer 
problem  has  been  addressed  by  a  variety  of  individuals.  Several  authors' *2,3  have  considered 
continuous  thrust  transfers  requiring  one  revolution  or  less  to  complete,  while  others4  have 
investigated  transfers  requiring  many  revolutions  using  some  version  of  the  method  of 
averaging.  This  lottar  approach  usually  leads  to  a  series  of  small,  nearly  impulsive,  burns  and 
tssuboptlmal. 

It  is  always  possible,  in  principle,  to  solve  the  transfer  problem  through  integration  of 
the  equations  of  motion.  Unfortunately,  the  amount  of  computation  required  for  transfers 
involving  many  revolutions  is  prohibitive,  and  numerical  errors  rapidly  become  unacceptably 
targe.  Wfesel  and  Alfaro?*6  have  developed  an  analytical  approach  which  leads  to  a  more 
tractable  solution  of  the  optimal  transfer  problem  for  many  revolutions.  Their  approach 
involves  splitting  the  control  problem  into  a  "fast"  timescale  problem  over  one  or  a  few  orbits, 
and  e  "slow"  timescale  problem  over  an  entire  transfer.  The  slow  timescale  solution 
incorporates  a  series  of  solutions  to  the  fast  timescale  problem.  This  approach  leads  to  an 
optimal  solution  for  an  infinite  number  of  revolutions.  Wiesel  and  Alfaro  found  that  the  errors 
in  obtaining  desired  final  conditions  were  proportional  to  I  /N,  where  N  was  the  number  of 
revolutions  in  the  transfer.  This  result  suggests  that  the  solution  of  the  optimal  control 
problem  for  finite  revolution  transfers  should  be  amenable  to  a  perturbation  theory  approach. 

This  paper  will  develop  an  analytic  solution  to  the  fast  timescale  problem  using  a 
first-order  perturbation  approach.  This  analytic  solution  is  compared  with  the  solution 
obtained  by  integrating  the  aquations  of  motion  to  determine  its  range  of  validity.  Baaed  on  these 
results,  some  conclusions  are  drawn  regarding  the  application  of  perturbation  theory  to  the  low 
thrust  transfer  problem. 


11.  Low  Thrust  Spiral  Solution 


Asa  limiting  case,  the  optimal  control  solution  for  a  vehicle  having  infinitesimally  small 
acceleration  is  developed  for  a  planar  transfer  beginning  in  e  circular  orbit.  The  solution  is 
optimal  In  the  sense  of  minimizing  the  time  required  to  perform  the  transfer  and,  hence, 
maximizing  the  time  rate  of  change  of  the  semi  major  axis,  a.  The  rate  of  change  of  the 
semimajor  axis  is  given  by7 


M  . 

dt 


0) 


where  the  energy  per  unit  mass  is 


(2) 


and  the  time  rate  of  change  of  the  energy  is 


d£ 

dt 


(3) 


In  the  above  aquations,  the  gravitational  parameter  is  represented  by  Ji,  m  is  the  vehicle 
mass,  v  is  the  true  anomaly,  and  F  and  Y  are  the  vehicle  thrust  and  velocity  vectors, 
respectively.  The  orbital  radius,  r,  is  the  same  as  the  semimajor  axis  for  a  circular  orbit  R 
and  T  are,  respectively,  the  radial  and  tangential  components  of  the  vehicle  acceleration.  The 
derivative  of  r  with  respect  to  v  is 


(4) 

&L  .  .  miQF. 
dv  t  +  a  coa  v 


where  a  is  the  orbital  eccentricity. 
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For  a  circular  orbit  eccentricity  is  zero,  so  the  right  hand  side  of  Equation  ( 4)  is  zero. 
Hence,  the  radial  term  in  Equation  (3)  is  initially  zero  and  the  vehicle  directs  its  thrust 
tangentially.  Since  the  vehicle  has  infinitesimally  small  acceleration,  orbital  eccentricity  must 
remain  infinitesimally  small  and  the  right  hand  side  of  Equation  ( 4)  will  always  be  nearly  zero. 
As  a  consequence,  in  accordance  with  Equation  (3),  a  radial  acceleration  component  will  not 
result  in  any  change  in  energy.  Since  energy  is  affected  only  by  the  tangential  acceleration 
component,  the  optimal  control  solution  is  to  direct  the  thrust  vector  tangent  to  the  orbit.  This 
corresponds  to  a  yaw  angle  of  zero. 

For  zero  yaw  angle,  the  dot  product  of  the  vehicle  thrust  vector  end  the  vehicle  velocity 
vector  is 


(S) 


FV-FV 


so  the  rate  of  change  of  the  semimejor  axis  is  given  by 


(6) 


Jfe  »  iflil v 

dt  H  m  ore  m 


where  vdrc  is  the  circular  orbit  velocity.  Separating  variables  and  integrating: 


(7) 
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yields  the  semimejor  axis  as  a  function  of  time: 


(8) 
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ii  VI  -  e2 


I  ejn  1/ 1 

fl  4  1 

I  dill  V  \ 

^  1  ♦  ecosv  '* 

jfi  _  r  sinu 


dt  na2Vl  -  e2  sin  i 


Unfortunately,  this  form  involves  the  argument  of  perigee,  <•>.  and  the  eccentricity,  e,  so  it  is 
singular  for  small  eccentricity.  This  singularity  problem  can  be  eliminated  by  making  the 
change  of  variables: 

(10) 

h  »  e  cosu 
k  =  e  sin  ia 


Another  problem  is  that  the  true  anomaly,  v,  is  not  well  defined  for  nearly  circular  orbits.  The 
true  anomaly  can  be  reptaoed  by  the  argument  of  latitude,  u,  where: 

(11) 

u  *  b>  ♦  v 

This  change  of  variables  transforms  the  independent  variable  from  time  ( implicit  in  v)  to  u. 


Applying  these  variable  changes  to  the  Lagrange  planetary  equations,  expanding  for  small 
eccentricity,  and  retaining  only  the  lowest  order  terms  yields:9 

02) 


N 


ilnw- 
n2a  sin  i 


N 


In  these  equations  R,  T,  and  N  are  the  acceleration  components  in  the  radial,  tangential,  and 
normal  directions,  respectively,  i  is  the  inclination,  and  Q  Is  the  right  ascension  of  the 
ascending  node.  The  quantity  a  is  the  mean  motion. 


(13) 

a"v/$ 

The  mass  and  acceleration  of  a  low  thrust  vehicle  making  small  orbit  changes  can  be 
treated  as  constant  This  assumption  can  not  be  made  when  considering  the  long  timescale 
problem.  For  a  vehicle  with  constant  acceleration  A,  pitch  angle  0,  and  yaw  angle  t  (see 
Figure  1 ),  the  acceleration  components  are: 


(14) 

T  -  A  cos  8  cos  ♦ 

R  ■  -A  cos  0  sin  ♦ 

N  -  A  sin  8 


,  . 
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Figure  1.  Vehicle  Acceleration  Components 


Another  useful  relationship  Is 


€  - 


(15) 


where  €  Is  the  dimensionless  ratio  of  vehicle  acceleration  to  local  gravitational  acceleration. 
Substituting  the  expressions  for  the  acceleration  components  and  €  Into  the  equations  of  motion 


H  ■  el  -  sinu  cosO  sin*  +  2  cosu  cos6  cost] 

^  «  elcosu  cosGsint  ♦  2  sinu  cos 8  cost! 

€  cosu  sinO 

■39  “  sinu  8in0 
du  sin  i 


Simplification  to  Pi«nr  Cmb 

It  was  decided  to  develop  only  the  more  limited  ceae  of  orbital  changes  within  e  single 
plane.  This  is  equivalent  to  setting  the  pitch  angle  6  equal  to  zero  in  the  equations  above.  This 
assumption  results  in  the  following  set  of  three  equations-. 


(17) 


-  €f-$inu  sint  ♦  2 cosu  cost] 


•  elcosu  sint  ♦  2  sinu  cost] 


This  is  the  final  form  of  the  equations  of  motion  used  to  develop  the  analytical  solution  to  the 
optimal  low  thrust  transfer  problem. 


To  solve  the  tow  thrust  transfer  problem ,  ft  fs  necessary  to  determine  the  optima)  yaw 
angle  program  to  perform  an  orbital  transfer  from  a  given  set  of  initial  conditions  (o0,  h0,  k^) 
to  a  given  set  of  final  conditions  ( a,.  hr,  k,).  For  a  vehicle  under  constant  thrust,  an  optimal 
trajectory  is  equivalent  to  a  minimum  time  trajectory.9  if  the  vehicle's  orbital  eooentricity  is 
very  small  at  all  points  in  the  trajectory  then  time  is  simply  proportional  to  u.  Based  on  this 
assumption,  the  Hamiltonian  for  the  optimal  transfer  is 

(18) 

Since  only  small  orbital  changes  are  being  considered,  the  dependant  variables  on  the  right  sides 
of  the  equations  of  motion  can  be  treated  as  constants.  Using  this  assumption,  the  canonical 
equations  for  the  lagrange  multipliers  are 

(19) 

-ft -o 

TT1  »  -#  -0 

ou  dk 

so  all  three  of  the  lagrange  multipliers  are  constants.  The  solution  to  the  long  timescale  problem 
would  incorporate  a  series  of  such  small  changes  Inserting  the  equations  of  motion  into  the 
Hamiltonian  and  computing  the  yaw  optimality  condition 


yields  the  control  law: 


tent  - 


-X2sinu  ♦  X3oosu 
2  ( a  X|  ♦  X2  cos  u  ♦  X3  sin  u ) 


Note  that,  even  though  the  X'sare  constant,  if  will  vary  because  of  its  dependence  on  u  The 
direction  of  the  transfer ,  inward  or  outward,  is  used  to  resolve  the  ambiguity  in  the  solution  for 
t  For  outward  transfers  ♦  will  take  on  values  In  the  range  -n/2  Or  i  if/2,  while  for 
inward  transfers  t  will  take  on  values  in  the  other  two  quadrants. 

Since  the  lagronge  multipliers  are  constant  throughout  a  given  orbital  transfer  so  long  as 
the  change  in  a  is  small,  and  ♦  is  dependent  on  these  multipliers,  the  multipliers  are 
effectively  the  control  variables  for  this  problem.  The  yaw  control  law  can  be  simplified  by 
letting 


U|"K? 


The  yaw  control  law  then  becomes. 


-UiSlnu  +  U2cosu 
2(q  +  Utcosu  +  U2s1nu ) 


To  further  simplify  the  analytical  solution  to  the  low  thrust  transfer  problem,  it  is 
assumed  diet  the  yaw  angle  t  is  small.  Since  the  optimal  yew  angle  for  a  vehicle  having 
infinitesimal  thrust  is  **0 ,  this  is  a  reasonable  assumption  for  a  low  thrust  vehicle.  Note  that 
if  ♦  is  smell,  U|  and  U2  must  be  small.  Using  this  small  angle  assumption,  the  equations  of 
motion  become 


“  €l-*sinu  ♦  2cosol 


&  “  €| tcosu  ♦  2sinu] 


and  the  yaw  control  law  becomes 

(25) 

.  -Uisinu  ♦  U2C0SU 

Y~  2a 

Note  that  the  equation  for  the  semimejor  axis  does  not  involve  ♦  At  this  level  of  appraximation 
the  analytical  solution  is  identical  to  the  spiral  solution  for  a.  The  differential  equation  tor  a 
can  be  solved  directly  in  terms  of  u  to  yield 


Q  -  Qf 


M  -  4Aa2(u  -u0) 


The  changes  in  orbital  elements  over  an  entire  orbit  transfer  are  computed  by  integrating 
the  equations  of  motion  over  u.  For  the  semimajor  axis  this  is 


AO.Of-Qo-j'-jjj 


and  similarly  for  the  other  two  orbital  elements.  Carrying  out  these  integrations  for  the 
equations  of  motion  resulting  from  the  small  ♦  assumption  yields 


\  £ 
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fi 
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(28) 


Ao  *  2do€u 


Ah  -  ^  [doSinu  ♦  ^-(-sinu  cosu  ♦  u)  -  •^■sin2u  | 

Ah  -  ■{£  (Oo(l  -  cos u )  ♦  ^sin2!*  ♦  ^-(sinu  cosu  ♦  u)| 


Note  that  Equation  (26)  is  valid  for  large  changes  in  a  while  the  first  of  Equations  (26)  is 
valid  only  for  small  changes. 


♦  ■lira. 


ol  Solution 

Given  the  initial  and  final  conditions  for  an  orbital  transfer  it  is  necessary  to  compute  the 
controls  ( i.e.  the  X's)  required  to  achieve  that  transfer.  It  is  assumed  that  each  transfer  begins 
and  ends  In  a  circular  orbit.  The  total  changes  in  h  and  k  over  an  entire  transfer  of  this  type 
must  therefore  be  zero,  since  these  two  elements  ere  directly  proportional  to  the  eccentricity. 
Using  this  fact,  the  variables  U,  and  ^  can  be  determined  by  simultaneously  solving  the 
following  two  aquations. 


(29) 

Ah  -  0  •  QoSinuf  ♦  -^(-sinufcosuf  *  Uf)  -  stn2Uf 
Ak  -  0  •  Op(l-  cosur)  ♦  ^sin2uf  ♦  sinuf  cosur  ♦  uf) 


Solution  of  these  equations  yields: 


U  i 


(30) 


(U2Sinuf  -  8q8)sinuf 
or  -  sinuf  cosuf 


eq.UloV  -  (1  -  oosuf)(ur  -  sinuf  cosuf)) 
°2  u?  -  sin2Uf  coe2Uf  ♦  sin*Uf 
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The  tegrange  multiplier  associated  with  the  semlmajor  axis  Is  computed  using  the 
traravarselity  condition 


H(.f)  ■  o  -  l  .  *u'&  .  u2£)| 


(31) 


Solving  this  equation  for  X«  yields 


(32) 


-2U 


*  Aa(4a(a  +  oosu  ♦  sinu)  ♦  (uisinu  -  U2oosu)(sinu  -  cosu)l 


where  a  as  a  function  of  u  is  given  by  Equation  26.  The  values  of  Xj  end  X3  are  cmputed 
from  X^UpandUj.  With  the  X's  in  hend,  the  yew  angle,  ♦.can  be  computed  for  any  value  of 
u  using  Equation  (25). 


,‘JI 


The  analytical  solution  to  the  low  thrust  transfer  problem  developed  in  Section  IV  was 
verified  by  comparison  with  the  integrated  equations  of  motion  (Eqns.  13)  for  a  set  of  orbital 
transfers.  All  computations  were  carried  out  in  canonical  units,  for  which  the  distance  unit  is 
one  Earth  radius  and  the  gravitational  parameter,  jl,  is  numerically  equal  to  one.  All  of  the 
verification  trajectories  began  at  u  -  0.  The  initial  and  final  conditions  for  all  of  these 
trajectories  were 

(33) 

(vW^t.o.o) 

(Ophpkf) -(1.01,0,0) 


Note  the  small  change  in  the  semimajor  axis.  Solutions  for  these  conditions  can  be  extended  to 
other  values  of  a,,  and  Of- by  scaling  £ 

Using  the  above  conditions  and  specifying  Uf,  the  X's  and  the  required  acceleration  were 
computed  In  order  to  make  a  valid  comparison,  these  computed  values  were  used  in  both  the 
analytical  solution  and  the  solution  via  direct  integration  of  the  aquations  of  motion.  For  the 
analytical  solution,  Equation  (26)  was  used  instead  of  the  first  of  Equations  (28)  to  provide 
higher  accuracy  for  changes  in  the  semimajor  axis.  In  the  integrated  solution,  however,  the 
analytical  result  of  Equation  ( 26)  was  not  used;  all  three  of  Equations  ( 1 7)  were  integrated 
Four  hummed  integration  steps  per  revolution  were  used  The  small  angle  assumption  was  not 
made  for  the  integrated  solution. 

Orbital  transfers  computed  using  both  the  analytical  solution  and  the  integrated  solution 
are  compared  in  Figures  2-9.  Figures  2-5  show  the  orbital  elements  and  yaw  angle  as  a 
function  of  u  for  a  transfer  to  ur  -  1.8  revolutions.  The  correspondence  between  the  analytic 
and  integrated  solutions  is  poor.  Note  that  for  small  changes  in  a  the  yaw  angle  computed  using 
the  two  solution  methods  is  essentially  identical ,  so  Figure  5  shows  only  one  curve.  Figures  6-9 
show  a  transfer  to  ur  -  49  revolutions,  and  in  this  case  the  correspondence  between  the  two 
solutions  is  very  good  Figures  2-9  indicate  that  the  analytics!  solution  is  valid  for  large  ur, 
but  breaks  down  at  small  values  of  uf  ( i.e.  higher  accelerations). 
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Using  the  analytical  solution,  the  maximum  yaw  angles  for  transfers  Involving  large 
values  of  iif-Cneer  100  revolutions)  were  computed.  Since  end  Of  are  fixed  (Eqns.  33), 
large  uf  corresponds  to  a  very  low  acceleration  level.  As  shown  In  Figure  10,  the  maximum 
yaw  angles  for  these  transfers  were  very  small.  Note  the  continuing  decrease  in  the  maximum 
value  of  the  yaw  angle  as  ur  Increases.  Hence,  for  very  tow  vehicle  accelerations  the  analytical 
solution  gives  results  close  to  those  of  the  spire)  solution. 
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Figure  8.  k  vs.  u  (uf  -  4.9  rev.) 
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A  solution  to  the  low  thrust  boundary  value  problem,  for  a  gtven  transfer,  occurs  when 
the  value  of  the  Hamiltonian  first  reaches  zero.  A  typical  plot  of  the  Hamiltonian  as  a  function  of 
u  is  shown  in  Figure  1 1.  The  Hamiltonian  oscilates  up  and  down,  but  its  mean  value  decreases 
os  the  semimajor  axis  increases.  If  the  Hamiltonian  is  extended  beyond  the  point  at  which  it 
first  equals  zero,  it  will  dip  below  zero  with  each  subsequent  revolution.  This  indicates  that  a 
minimum  fuel  solution  to  die  low  thrust  transfer  problem,  as  opposed  to  minimum  time 
solution,  can  be  developed  For  a  minimum  fuel  solution,  the  vehicle's  engine  would  be  rm  when 
the  Hamiltonian  was  positive  and  off  when  it  was  negative.  The  turn  on  and  turn  off  point*  for  an 
Inward  transfer  from  ur  to  u0  would  simply  be  the  mirror  image  of  the  turn  on  and  turn  off 
points  for  an  outward  transfer  from  u0  to  ur 

The  quantities  €,  Xt,  X3,  and  the  absolute  value  of  1'max  are  shown  as  a  function  of 
uf  in  Figures  12-16.  Note  that  solutions  do  not  exist  for  some  values  of  uf  below  2.5 
revolutions.  All  three  X's  are  periodic  with  a  period  of  2n  (one  revolution)  and.  in  each 
revolution,  they  are  all  singular  at  the  0.5  and  0.75  revolution  points.  Because  of  dependence 
on  the  inverse  rate  of  change  of  the  semimajor  axis,  the  mean  magnitude  of  X,  increases  as  uf 
increases.  The  mean  values  of  Xj  and  X3  do  not  change.  also  has  a  period  of  2n,but 
its  maximum  value  per  revolution  asymptotically  approaches  zero  as  ur  Increases  and  vehicle 
acceleration  decreases.  The  decrease  in  as  ur  increases  is  shown  in  both  Figure  16  and 
Figure  10. 

One  major  advantage  of  the  analytical  solution  over  the  integrated  solution  is  that  it 
requires  much  less  computation.  For  any  transfer,  using  the  analytical  solution,  the  values  of 
the  orbital  elements  for  any  given  value  of  u  can  be  computed  directly  from  Equations  ( 26).  In 
comparison,  obtaining  the  orbital  elements  for  the  same  value  of  u  using  the  integrated  solution 
requires  integrating  over  the  entire  range  from  u  ■  0  to  u.  Computing  and  outputting  the 
orbital  elements  forty  times  per  revolution,  the  analytical  solution  was  able  to  complete  one 
revolution  about  ninety  times  as  quickly  as  the  integrated  solution.  This  gives  a  rough  indication 
of  the  computational  advantage  of  the  analytical  solution. 


VII.  Conclusions  and  Reo 
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In  this  paper  an  analytic  solution  to  the  optimal  low  thrust  transfer  problem  was 
developed  using  a  first-ordBT  perturbation  approach.  For  low  vehicle  accelerations  this 
solution  yielded  results  quite  close  to  those  obtained  by  integrating  the  vehicle  equations  of 
motion.  At  higher  acceleration  levels,  however,  the  analytic  solution  diverged  from  the 
integrated  solution.  For  very  low  acceleration  levels  the  analytic  solution  was  found  to  approach 
the  limiting  case  of  the  infinitesimally  low  thrust  spiral  solution. 

The  results  described  in  this  paper  could  be  further  developed  by  incorporating  a  higher 
order  perturbation  model.  Such  a  higher  order  model  would  extend  the  validity  of  the  analytic 
solution  to  higher  levels  of  vehicle  acceleration.  Another  useful  activity  would  be  to  develop  an 
analytical  solution  for  the  full,  three  dimensional  orbital  transfer  problem,  eliminating  the 
limitation  to  planar  transfers.  Using  the  minimum  time  transfer  solution  as  a  basis,  the 
minimum  fuel  solution  should  also  be  developed.  By  incorporating  results  from  all  of  these 
recommended  areas  of  stud/,  there  is  potential  for  the  development  of  powerful  and  efficient 
tools  for  the  analysis  of  low  thrust  transfers. 
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